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In classical superconductors an energy gap and phase coherence appear simultaneously with pairing at the 
transition to the superconducting state. In high-temperature superconductors, the possibility that pairing 
and phase coherence are distinct and independent processes has led to intense experimental search of their 
separate manifestations. Using femtosecond spectroscopy methods we now show that it is possible to clearly 
separate fluctuation dynamics of the superconducting pairing amplitude from the phase relaxation above 
the critical transition temperature. Empirically establishing a close correspondence between the superfluid 
density measured by THz spectroscopy and superconducting optical pump-probe response over a wide 
region of temperature, we find that in differently doped Bi2Sr 2 CaCu 2 O s +j crystals the pairing gap 
amplitude monotonically extends well beyond T a while the phase coherence shows a pronounced power-law 
divergence as T— > T a thus showing that phase coherence and gap formation are distinct processes which 
occur on different timescales. 

Anomalous normal state behavior above the critical temperature appears to be a hallmark of unconven- 
tional superconductivity and is present in many different classes of materials. A pseudogap (PG) state has 
been suggested to be associated with a wide range of possible phenomena preceding the onset of mac- 
roscopic phase coherence at the superconducting (SC) critical transition temperature at T c : pre-formed pairs 1 " 9 , a 
spin-gap 10 , the formation of a Bose metal 11 , a Fermi or Bose glass, or a state composed of "dirty bosons" 12 " 14 , and 
more recently a charge-density-wave state 1516 . 

In addition to the PG response below the temperature designated as T*, the response attributed to "super- 
conducting fluctuations" above T c has been observed in a number of experiments 17 " 27 . The temperature region T c 
<T< T onset where such fluctuations are observable is significantly wider than in conventional superconductors, 
but smaller than T*. The open and obvious question is whether the pseudogap, or the superconducting fluctua- 
tions can be attributed to pairing. 

The problem in separating the response due to superconducting fluctuations from the PG is that so far, 
inevitably, one has had to make extrapolations, or assumptions about the response functions underlying tem- 
perature dependences and line shapes in transport 1718 ' 28 " 30 , magnetic susceptibility 26,31 , specific heat 20,21 or photo- 
emission (ARPES) 22 , which may at best introduce inaccuracies in the temperature scales, and at worst lead to 
erroneous conclusions. Alternatively one can suppress superconductivity by high magnetic fields up to 60 T 32 , 
although there exists a risk of inducing new states by such a high field 33 . Thus, so far it has not been possible to 
satisfactorily characterize superconducting fluctuations and discriminate between fluctuations of the amplitude 
d\j/ (related to the pairing gap) and phase 30 of the complex order parameter = \jje id . 

In pump-probe experiments three relaxation components shown in Fig. 1 a) are typically observed 34 : 1) the 
quasiparticle (QP) recombination in the SC state, 2) pseudogap state response below T* and 3) energy relaxation 
of hot electrons. The QP dynamics has been shown to be described very well by the Rothwarf-Taylor (R-T) 
model 35,36 , and the response related to the presence of non-equilibrium QPs is thus unambiguous. Importantly, 
the presence of the QP response is directly related to the presence of a pairing gap for QP excitations. 

Pump-probe experiments have already shown the coexistence of the pseudogap excitations with supercon- 
ductivity below T c over the entire range of phase diagram 37 " 41 . However, little attention has been paid to the region 
of superconducting fluctuation. In this paper we present measurements by a 3 -pulse technique which allows us to 
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Figure 1 | Description of the 3 pulse method and pump-probe signals. Underdoped sample data: (a) Transient reflectivity signals observed at different 
temperatures in Bi-2212 - SC signal, the PG signal and energy relaxation of hot carriers (Inset: susceptibility curve shows T c = 81 K for the underdoped 
sample), (b) The QP recombination at different T above T c . (c) Two typical 3-pulse pump-probe traces at f^pO.2 and 4.5 ps show signal with suppressed 
and recovered superconducting component respectively, (d) Pulse sequence and delays notation in the 3-pulse experiment, (e) The pump-probe response 
at 120 K in the absence of the D-pulse and 0.2 ps after the D pulse arrival. (The noise level is larger in comparison to c) due to a shorter signal 
accumulation time.) (f) A schematic representation of the 3-pulse experiment. Colors are schematic. Destruction (D) pulse destroys the 
superconductivity, Pump-probe sequence probes the recovery of the quasiparticle response, (g) A typical result of three pulse experiments at 90K shows 
suppression of the SC component after the D pulse arrival and gradual recovery with tp^p. Readings along the blue and dark- red line are shown in Fig. lc). 
The color represents amplitude of the reflectivity change. The values of the color bars indicate AR/R X 10~ 5 . (h) Recovery of the superconducting 
component amplitude with tpj_p. 



single out the response of superconducting gap fluctuations. The 
approach is based on selective destruction of the superconducting 
state by a femtosecond laser pulse 42 which allows us to discriminate 
pseudogap excitations from superconducting fluctuation seen in 
transient reflectivity signals, thus avoiding the necessity of making 
extrapolations or assumptions in separating the different contribu- 
tions. We then compare these data with a.c. conductivity (THz) 
measurements from 43 and establish proportionality of the amplitude 
of superconducting component in pump-probe experiment to the 
bare phase stiffness p 0 > measured by THz experiments. This is 
directly proportional to the bare pair density /? 5 23,43 , which in turn 
coincides with the superfluid density when the latter is measured on a 
timescale on which changes of the order parameter due to either de- 
pairing or movement of the vortices can be neglected. Comparison of 
the critical behavior of the amplitude and phase correlation times 
above T c leads us to the conclusion that two quantities arise from 
different microscopic processes. 

We perform measurements on under- (UD), near optimally- (OP) 
and over-(OD) doped Bi 2 Sr 2 CaCu 2 0 8+ s (Bi2212) with T c s of 81, 85 
and 80 K respectively. In the discussion we focus on the underdoped 
sample, and discuss comparisons with the optimally and overdoped 
samples, where applicable. 

Results 

Measurements of pairing amplitude above T c . To separate the SC 
component from the PG component we use a 3 pulse technique 
described in Refs. 44-47, and schematically represented in Fig. If). 



A pulse train of 800 nm 50 fs pulses produced by a 250 kHz 
regenerative amplifier is divided into three beams with variable 
delays. First a relatively strong "destruction" (D) pulse, with fluence 
just above the superconducting state photodestruction threshold, 
F^= 13 fi] I 'cm 2 48 , destroys the superconducting condensate 49 . The 
ensuing recovery of the signal is measured by means of the 2 pulse 
Pump-probe (P-pr) response at a variable delay t D - P between D and 
P pulses. There is no effect of the D pulse on the PG response within 
the noise level of <jar/r <0.4x10 -5 indicated by measurements at 
120 K shown in Fig. le). A typical result of the 3 pulse experiment is 
presented in Fig. lg). In the absence of the D pulse the signal consists 
of a positive SC and a negative PG component. After the arrival of the 
D pulse we see a disappearance of the SC part and only the PG 
component is present (dark- red line in Fig. lc)). With increasing 
t DP the superconducting response gradually re- emerges (blue line 
on Fig. lc)). As most of the condensate in the probe volume is "des- 
troyed" by the D pulse we can extract the superconducting 
component by subtracting the signal remanent after the destruction 
(measured 200 fs after the D pulse) from the signal obtained in the 
absence of the D pulse. Such an extracted superconducting 
component is plotted in Fig. lb). The signal is detectable up to 
T onset = 104 K, which is 0.28 T c above T c but much lower than 
T* ~ 2.5 T c . In Fig. lh) we show the recovery of the amplitude of 
the superconducting component A sc as a function of the time delay 
t D - P after the D pulse for 90 K. The temperature dependence of the 
amplitude of the SC component measured by 3-pulse technique 

A 3pulse is shown m pig 2b) 
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Figure 2 | Comparison of pairing amplitude and phase coherence 
dynamics, (a) The recovery time of the optical superconducting signal 
(T rec ), the QP recombination time z 3 Q P ulse measured by the three pulse 
technique and the QP recombination time T s ^ traction from two pulse 
Pump-probe measurements obtained by subtraction of the PG signal. A fit 
to the data using a BKT model (eq. 3 of ref. 19) is shown by the solid red 
line. The dashed line shows the fluctuation lifetime z GL given by the TDGL 
theory. The phase correlation time z™ z obtained from the THz 
conductivity measurements 43 is also shown for comparison, (b) The 
amplitude of the SC signal measured by the three pulse technique (A 3 s ^ ulse ) 
and the two pulse measurements with the PG signal subtracted 

^Subtraction^ ^ ^ phase stiffness p( 43 ( norma l ize d at Tc ) shows 

remarkable agreement with the optical response. 



For comparison, standard pump-probe measurements need to 
separate the SC relaxation from the PG relaxation by subtraction 
of the high temperature response extrapolated into the supercon- 
ducting region. This approach suffers from the same uncertainties 
as other techniques such as conductivity, heat capacity, diamagnet- 
ism and ARPES. The actual T-dependence of the PG response can 
vary with doping, pump fluence and probe wavelength 48,50 . But, in 
Fig. 2b) we show that the subtraction procedure - with the use of a 
model 48 - gives results in agreement with the direct 3 pulse measure- 
ments. The remaining discrepancies in the amplitude can be 



explained by an incomplete destruction of fluctuating superconduct- 
ing state in the 3 pulse experiment and errors in the PG subtraction. 

In Fig. 2a) we see that the QP relaxation time t 3 q P u1s6 obtained by 
fitting an exponential function to the data Fig. lb) decreases rather 
gradually with increasing T above T c and nearly coincide with the QP 
relaxation time obtained by the pseudogap subtraction procedure 
Subtraction _ recovery time x rec obtained from exponential fits to 
the recovery of the SC response above T c (Fig. 1 h) shows a similar In- 
dependence. The experiments thus show that the recovery of the SC 
gap and the QP relaxation show very similar dynamics above T c . 

Comparison of optical and a.c. conductivity measurements. We 

now compare these data with THz measurements of the order 
parameter correlation time and bare phase stiffness 19 ' 43 . The 
agreement between p 0 and A sc shown in Fig. 2b) is seen to be 
remarkably good over the entire range of measurements 0.8 T c < 
T < 1.3 T c . This agreement is important because, taking into account 
n s ~ |*F| 2 , it validates the approximation that the pump -induced 
changes in the reflectivity or dielectric constant e for small n s are 
related to the order parameter *F as 5R ~ ds ~ I^F] 2 . 

In contrast to A sc and p 0 , remarkable differences are seen in the 
temperature dependences of the characteristic lifetimes shown in 
Fig. 2a) obtained by optical techniques and THz conductivity mea- 
surements. The phase correlation time Tq Hz determined from the 
THz conductivity 19 dies out very rapidly with increasing temper- 
ature, while the T-dependence of the amplitude relaxation (T 3 Qp lse , 

Subtraction Qr {$ much mQre gmduaL 

Measurements on an optimally doped sample (Fig. 3c, d))) show 
qualitatively the same results with T onset ~ 102 K, which is 17 K 
above T c , and slightly faster decrease of both amplitude and QP 
relaxation time with temperature. For the overdoped sample, F^ G 
becomes comparable to F^ c , so the superconducting component 
cannot be significantly suppressed without affecting the pseudogap. 
Nevertheless the superconducting component is clearly observable in 
the 2-pulse response up to T onset ~ 93 K, i.e. 13 K above T c . 
Comparison of 2-pulse data for different doping levels is shown in 
Fig. 3 e)-g), showing qualitatively similar behavior of the SC ampli- 
tude above T r . 



Discussion 

The co-existence and distinct dynamics of the PG and SC excitations 
above T c highlights the highly unconventional nature of these states. 
A possible explanation for the coexistence of the SC and PG excita- 
tions is that the SC and PG quasiparticles which are giving rise to the 
observed processes are associated with relaxation at different regions 
on the Fermi surface. Recent Raman and cellular dynamical mean- 



3 

IF o 
o. 2. 

£ CL A 
00 O ' 

0 
10 



q- o 

< 



a)UD 

K - 


c) OP 
■ ■ 


b)UD 
■ 


_d)OP 
■■ 




0 5 10 15 20 0 5 10 15 20 
T-T C [K] T-T C [K] 

Figure 3 | Doping dependence. Comparison of the QP recombination time T 3 Qp lse (a) and (c) and the amplitude of the SC signal (b) and (d) 
measured by the three pulse technique for under-(UD) and optimally(OP) doped samples, respectively, (e)-(g) T-dependence of 2-pulse response for 
under- (UD), optimally (OP) and overdoped(OD) samples. Blue dashed and cyan solid lines marks T onset and T c respectively. T onset shows gradual 
decrease with doping. The values of the color bars indicate ARIR X 10~ 4 . 
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field studies 51 have suggested that the PG may originate from the 
states inside the Mott gap, which are characterized by s-wave 
symmetry and very weak dispersion. Such a localized nature of 
the PG state excitations is consistent with previous assignments 
made on the basis of pump-probe experiments 39,48 . In contrast, the 
superconducting gap fluctuations have predominantly d-wave 
symmetry 52 and are more delocalized. This would explain the sim- 
ultaneous presence of the SC fluctuations and PG components in 
pump-probe experiments. 

Perhaps the most widely discussed model in the context of distinct 
pairing and phase coherence phenomena is the Berezinskii- 
Kosterlitz-Thouless (BKT) transition 53 " 55 by which decreasing tem- 
perature through T onset and approaching T c causes freely moving 
thermally activated vortices and anti-vortices to form pairs, thus 
allowing the condensate to acquire long range phase coherence in 
an infinite-order phase transition sharply at T c . The bare pair density 
is finite to much higher temperatures (up to T onset ), where pairing is 
caused by a different mechanism, and the pseudogap is considered to 
be an unrelated phenomenon 56 " 59 . The effect is characterised in terms 
of a phase stiffness p s , a quantity which characterizes the destruction 
of phase coherence by thermal fluctuations at a temperature T c = 
Tbkt — np s /8. In cuprates p s is small due to reduced dimensionality 
and the low carrier density. Within this approach a.c. (THz) con- 
ductivity 19 ' 43 , heat capacity 20,21 , diamagnetism 26,31 , ARPES 22 and 
Nernst effect 23 " 25,60 measurements were interpreted. 

The BSCCO family is the most two-dimensional of the cuprate 
materials, so here the BKT mechanism for describing the order above 
T c would be expected to be most applicable 43 . However, the data in 
Fig. 2a) show that the drop in z e is not nearly as abrupt as the BKT 
model predicts. One possibility is that this broadening arises from 
chemical inhomogeneity of the sample, but the absence of a peak in 
the heat capacity at T omet 20 also appears to exclude the possibility of a 
pure BKT transition, and implies amplitude fluctuations might be 
present between T c and T onset as well 21 . Thus additional mechanisms 
beyond BKT may also be present which would broaden the phase 
coherence transition, such as interlayer phase fluctuations. In this 
case the observed T-dependence would reflect the interlayer de- 
coherence 19,43 . 

In more traditional approaches using time- dependent Ginzburg- 
Landau (TDGL) theory 61 , thermal fluctuations are small for tempera- 
tures higher than ~2 K above T c 62 , but can give an observable con- 
tribution to the conductivity in this temperature range 17,18,28 " 30 . 
Relaxation within TDGL theory has a "longitudinal" relaxation time 
Ta, which corresponds to the relaxation of the magnitude of the SC 
order parameter, and a "transverse" relaxation time tq which corre- 
sponds to the relaxation of its phase. They are related to each other in 
magnitude, but have the same temperature dependence near T c , 
namely t gl ~ 1/(T — T c ) 63 . Perhaps unexpectedly, the temperature 
dependence of t™ z nearly coincides with the behavior predicted for 
tq by time- dependent Ginzburg Landau theory for Gaussian fluctua- 
tions. Fig. 2a) suggests that phase coherence within this system is 
established in a narrow, ~5 K temperature interval. However, the 
distinctly different critical behavior of the pair amplitude dynamics 
speaks in favor of unconventional models of superconductivity in 
which pairing and phase coherence occur independently, by different 
mechanisms. The implication is that the observed pairing amplitude 
which extends to more than 25 K above T c reflects the response of an 
inhomogeneous ensemble of gapped patches which are not mutually 
phase coherent. The weak temperature dependence of the amplitude 
cannot be described either by TDGL or BKT models. 

Beyond the BKT vortex and TDGL scenario, other phase-locking 
scenarios, such as Bose-Einstein condensation of bipolarons 2,3,64 and 
phase-coherence percolation 4,5 may also be consistent with the 
observed behaviour. In both of these cases pairing and phase coher- 
ence are also distinct processes. The former comes from the con- 
densation of pairs at T c as preformed pair kinetic energy is reduced, 



while percolation dynamics is associated with the time dynamics of 
Josephson tunneling between fluctuating pairs or superconducting 
patches. The percolation timescale t 7 is given by the Josephson 
energy Ej = I c ^ 0 /27r, where I c is the critical current and 0 O is tne flux 
quantum. In cuprates, tj = H/Ej~300 fs, which is compatible with 
the dynamics of phase shown in Fig. 2a). 

A picture highlighted by Fig. 2 thus emerges in these materials 
where the relaxation of the phase 0 is faster than relaxation of the 
amplitude ijj of the complex order parameter *F = \j/e i6 , the dynamics 
of \jj and 0 being governed by microscopically different processes. It is 
worth remarking here that the opposite situation is found in charge 
density wave dynamics, where the phase relaxation is slow compared 
to the amplitude relaxation, and the dynamics can be described by 
TDGL equations for the amplitude neglecting phase relaxation 0 45 . 

Methods 

Samples. The samples used in this work were under-, near optimally- and over- 
doped Bi2212 with T c s of 81, 85 and 80 K respectively. Samples were grown by the 
traveling solvent floating zone method. Critical temperatures were obtained from 
susceptibility measurements (e.g. inset in Fig. la) for the underdoped sample). 
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